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Unusual Quasiparticle Tunneling in High- Tc Cuprate Superconductors: Evidence for 
the BCS and Polaronic Multi-Gap Effects on Tunneling Spectra 

S. Dzhumanov0 O.K. Ganiev, and Sh.S. Djumanov 

Institute of Nuclear Physics, Uzbek Academy of Sciences, 100214, Ulughbek, Tashkent, Uzbekistan 

We propose a model of quasiparticle tunneling across the high- Tc superconductor-insulator-normal 
metal junction based on the different mechanisms for tunneling of electrons at positive bias and 
dissociating polaronic Cooper pairs and large polarons at negative bias, and the gap inhomogeneity 
(i.e., multi-gap) picture. We show that the main features of the tunneling spectra such as low- 
bias U- and V-shaped features, asymmetry and high-bias dip-hump features, their temperature and 
doping dependences, and shoulders inside the conductance peaks observed in high-Tc cuprates arise 
naturally from the model. The experimental tunneling spectra of Ba2Sr2CaCu2 08+« are fitted quite 
well by taking into account the distribution of BCS and polaronic gap values. 

PACS numbers: 71.38.+i, 74.20. Fg, 74.50.-|-r, 74.72.-h 



In conventional superconductors, superconductivity 
arises from the binding of electrons into Cooper pairs and 
the BCS gap opening in the electronic excitation spec- 
trum at the superconducting (SC) transition temperature 
Tc serves as the SC order parameter. In contrast, in high- 
Tc cuprate superconductors the precursor Cooper pairing 
of carriers occurs in the normal state and the BCS-like 
gap is manifested as a pseudogap (PC) opening at a tem- 
perature T* higher than Tc at which the Cooper pairs 
condense into a superfluid (SF) Bose-liquid state [H, 
As argued in Refs. UHl], in these materials the Cooper 
pairing is only a necessary but not a sufficient condition 
for the occurrence of superconductivity and the BCS-like 
gap might be different from the SC order parameter ap- 
pearing below Tc- Other alternative PG scenarios were 
also proposed (for a review, see, e.g., Refs. [341|)> among 
which the SC fluctuation scenario suggests that the PG 
appearing below T* is related to superconductivity as a 
form of precursor pairing and evolves smoothly into the 
SC gap at Tc Q. It is still highly debated whether the 
PG observed by various experimental techniques origi- 
nates from SC fluctuations above Tc or from the 
precursor non-SC Cooper pai ring below T* or from 
some other effects (see Refs. [iMl)- 

Scanning tunneling microscopy and spectroscopy 
(STM and STS) ^^U^ and angle-resolved photoemis- 
sion spectroscopy (ARPES) |18l - [2l| have made signifi- 
cant progress in the studies of PG phenomena in high-Tc 
cuprates and other materials. Progressive investigations 
of tunneling and ARPES spectra have provided impor- 
tant information on single-particle excitation gaps in the 
cuprates. The STM and STS techniques are very sen- 
sitive to the quasiparticle density of states (DOS), with 
the unique capability to measure any excitation gaps at 
the Fermi energy, and to the electronic (or gap) inhomo- 
geneities that are intrinsic to the cuprates [iB - flTl [23 - |2^ . 
These tunneling studies of the cuprate superconductors 
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have revealed a rich variety of tunneling spectra (differ- 
ential conductance- voltage {dl/dV — V) characteristics) 
of high- Tc cuprate superconductor (IITSC)-insulator (I)- 
normal metal (N) (called SIN) junctions. Two distinc- 
tive features of the dl/dV spectra are nearly U- and V- 
shaped characteristics observed in tunneling experiments 
on high- Tc cuprates [l^ [23 - l29j . Other distinctive tun- 
neling features systematically observed in the dl/dV — V 
characteristics of SIN junctions are:(i) asymmetric con- 
ductance peaks [26l - l30j , (ii) dip-hump structures appear- 
ing outside the conductance peak on the negative bias 
side [H, [13, m, [13], (iii) suppression of the peak on 
the negative bias side with increasing temperature and 
its vanishing somewhat below Tc or near Tc, leaving the 
hump feature (i.e., linearly increasing conductance) and 
the second conductance peak (on the positive bias side) 
O [271 [^. and (iv) shoulders inside the conductance 
peaks |22l424| . Similar peak-dip- hump feature and its 
persistence above Tc were also observed in ARPES spec- 
tra [11 [H. 

The unusual features of the tunneling spectra of high- 
Tc cuprates are neither expected within the simple s-wave 
BCS model nor within the BCS models based on the d- 
wave gap symmetry. For example, the tunneling spectra 
showing more flatter, BCS-like feature is difficult to re- 
produce in the d-wave model. While the other type of 
tunneling spectra showing more V-shaped feature might 
be expected either within the d-wave gap model (which 
fails to reproduce quantitatively the conductance peak 
height and shape [H [H, [l^l ) or within the s-wave multi- 
gap model proposed in the present work. The origins 
of the peak-dip-hump feature and the asymmetry of the 
conductance peaks have also been the subject of con- 
troversy. These features of the tunneling spectra have 
been attributed either to extrinsic (band-structure) ef- 
fects (i.e., the van Hove singularity and bilayer splitting) 
(3^ [33I or to intrinsic effects such as particle-hole asym- 
metry [13, [11] and strong coupling effects (which origi- 
nate from the coupling to a phonon mode j3^] or to a col- 
lective electronic mode [13] )■ Although several theoreti- 
cal models were successful in reproducing some tunneling 
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spectra with the peak-dip-hump features and asymmetric 
or nearly symmetric peaks observed in cuprates, the weh- 
established experimental tunneling spectra with different 
peak-dip-hump features [H, Ull (e.g., high-bias conduc- 
tances, which are nearly flat, linearly increasing at nega- 
tive bias or decreasing at positive bias, temperature- and 
doping-dependent peaks, dip-hump features and asym- 
metry of the conductance peaks) characteristic of high- 
Tc cuprates were not explained yet. The quantitative fits 
of the experimental spectra of Ba2Sr2CaCu20g+5 (Bi- 
2212) to a BCS d-wave gap and a T-independent gap 
functions IS] have been unsuccessful in reproducing the 
asymmetry (which is opposite to that of the observed 
tunneling spectra) and the temperature-dependent con- 
ductance curves. For underdoped Bi-2212, the asymme- 
try of the conductance peaks and its doping dependence 
found in Ref. E3(| are also inconsistent (i.e., in contrast) 
with other experimental data [H, [13, HI] . Further, STM 
and STS studies have shown that the electronic system 
of Bi-2212 is inhomogeneous and the gap distribution 
(i.e., inhomogeneity) has a strong effect on the tunneling 
spectra [22 24 Hj- It is necessary for any theoretical 
model to explain not only the asymmetry and peak-dip- 
hump features in tunneling conductance but also their 
evolution with temperature and doping, a flat (i.e., Li- 
shaped) and more V- shaped subgap conductances, the 
gap-inhomogeneity-induced shoulders inside the conduc- 
tance peaks that are observed recently by STM in tun- 
neling spectra of Bi-2212 [ll,!!!]. 

In this paper, we propose a simple and quite effec- 
tive model of quasiparticle tunneling based on the dif- 
ferent mechanisms for tunneling of charge carriers across 
the SIN junction at negative and positive biases and the 
multi-gap (i.e., gap inhomogeneity) picture. The pro- 
posed model reproduces the well-established experimen- 
tal tunneling spectra of high- Tc cuprates and their nearly 
U and V-shaped features, asymmetry, peak-dip-hump 
structure and shoulder-like features inside the conduc- 
tance peaks. We focus on the Bi-2212 system (which 
has been well studied experimentally) and show that the 
main experimental features of the tunneling spectra and 
their temperature and doping dependences can be well 
reproduced by using a BCS DOS at positive bias voltages 
{V > 0) and the combined BCS DOS and quasi-free-state 
DOS (originating from the dissociation of large polarons) 
at negative bias voltages {V < 0), and taking into ac- 
count the distribution of BCS and polaronic gap values. 
There is now ample reason to believe that the electron- 
phonon interaction in cuprates is strong enough and the 
relevant charge carriers in these systems are large po- 
larons [3, [S^ll^. Therefore, the precursor Cooper pair- 
ing of large polarons may occur above Tc with opening 
the energy gap A in their excitation spectrum [H-Q. As 
argued in Ref. [i^l , the binding energies of large polarons 
Ap and Cooper-like large polaron pairs are manifested as 
the two distinct non-SC gaps in high-T^ cuprates, one a 
temperature independent PC and the other a BCS-like 
gap. Two gap-like features observed in ARPES and tun- 



neling experiments are often misinterpreted as the coex- 
isting PC and SC gap. 

Unconventional SIN tunneling. - We consider the 
model which describes two specific mechanisms for quasi- 
particle tunneling across the SIN junction at < and 
V > 0, and explains the asymmetry of the tunneling 
current taking into account the different tunneling DOS 
existing in these cases. The first mechanism describes 
the S ^ N tunneling processes associated with the dis- 
sociation of Cooper-like polaron pairs and large polarons 
at < 0. In this case the Cooper pair dissociates into 
an electron in a normal metal and a polaron in a polaron 
band of the HTSC. This S N tunneling is allowed 
only at \eV\ > A. The dissociation of large polaron oc- 
curs at \eV\ > Ap and the carrier released from the po- 
laron potential well can tunnel from the quasi-free state 
into the free states of the normal metal. Such a. S —?' N 
transition gives an additional contribution to the tunnel- 
ing current. The other mechanism describes the electron 
tunneling from the normal metal to the BCS-like quasi- 
particle states in HTSC at V > 0, while the quasi- free 
states appearing only at the polaron dissociation are ab- 
sent. Therefore, at V > the tunneling current across 
SIN junction is proportional to the BCS DOS of HTSC 
given by 



DBcsiE,A) = 







for \E\ < A, 



where D{eF) = rrC'J'^ ^fep I ^f^'i^'^'fv' is the normal state 
DOS, rrip and sp are the mass and Fermi energy of large 
polarons. 

In the case < 0, the total current is the sum of 
two tunneling currents and proportional to the square of 
the tunneling matrix element, |Af p [4l|, the Dbcs{E, A) 
and the quasi-free state DOS. This current fiows from 
HTSC to normal metal at the dissociation of Cooper 
pairs and large polarons. In HTSC, the quasi-free carri- 
ers appearing at the dissociation of large polarons have 
the effective mass to* and energy E = Ap + H'^k'^ /2m*. 
Then the quasi-free state DOS is defined as 

Df(E AJ = / D{ei)^{\E\-Ap)/ei for \E\ > Ap, .3) 
1 for|i?| <Ap, 

where Sp is the Fermi energy of quasi-free carriers, 

D{ep) — (m*)^/'^ ^J^/V^-K^h^ . For the normal metal, 
the DOS at the Fermi energy Ep is independent of energy 
E, i.e., D{E) ~ D{Ef). Thus, at V > the tunneling 
current from the normal metal to HTSC is 

iN^siV) = C\M\^D{EF)D{eF) 

+ 00 

^{E + eV f - A2 

+ 00 

f /^J^[/(^-eV^)-/(e)]ds, (3) 



3 




100 
Voltage (mV) 



200 



FIG. 1: Main panel: SIN tunneling conductance for T=40 K calcu- 
lated using the one-gap model with single s-wave BCS gap A=35 
meV and single polaronic gap Ap=40 meV, exhibiting U-shaped 
feature at low-bias. Inset: comparison of the model (A=28 meV 
and Ap=22 meV, solid Une) with optimally doped Bi-2212 (Tc=92 
K) tunneling data at 4.8 K (dashed line) [lal . 



where G = eC\M\'^D{EF)D{eF), C is a constant, /(e) 
is the Fermi function, e = E + eV. The differential con- 
ductance, dI]sf-^s/dV is then given by 

dlN^s/dV = G{AiiAT,av)+A2iAT,av)), (4) 

where 

+00 

X exp[— a; — av]dx 



Ai{AT,av) 



A2{AT,av) 



+00 



^ — A|,(exp[— a; — ay\ + 1)^ ' 

X exp[a; — ay\dx 
\/ x^ — A|,(exp[x — av\ + 1)^ 



x — e/kBT, ay = eV/kBT. At negative bias voltages 
V < 0, the tunneling current and differential conductance 
are given by 

{+00 
— 00 



D{£F)y/e^ 



+ 00 



^\e\-Ap[f{e)-f{e + eV)]de\, (5) 



and 



dls^N 



dV 



= G{Ai{AT,'av) + A2{AT,-av) 
+aFiT)[Bi{A;,av) + B2iA;,av)]} , 

where e = E — eV, 

00 

ui\* \ f fTi a7 exp[x + av]dx 
B^{A^,ay) ^lj^^^^^^^Y^2^ 



(6) 



00 



S2(A;,ay) 



x\ - A 



exp[— a; + av]dx 
'P (exp[-x -t- ay] + ly 
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FIG. 2: Tunneling conductance as a function of temperature (a) 
and doping (b), calculated using the three-gap model. In (a) the 
gap values are: A=40, 30 and 22 meV, Ap=65, 47 and 31 meV. In 



(b) T=40 K and set of gap values: 
38 and 28 meV for curve 1: A=30, 



22 meV for curve 2; 
meV for curve 3. 



A=40, 30 and 24 meV, Ap=52, 
24 and 18 meV, Ap=33, 28 and 
and A=25, 20 and 16 meV, A„=21, 17 and 14 



aFiT) = [Diei)/D{sF)WkBT/eF, 



A* — A 



The SIN tunneling conductance curve calculated at r=30 
K for the single-gap case (concerning both the pola- 
ronic gap Ap and the s-wave BCS gap A) is shown in 
Fig. 1 (main panel). In this simple model, the absence 
of gap distribution would lead to the U-shaped spectral 
behavior at low bias and such a more flatter subgap con- 
ductance would be expected for homogeneous high- Tc 
cuprates. As can be seen in Fig.l, there are dip-hump 
feature and asymmetric peaks, with the higher peak in 
the negative bias voltage. The model is compared with 
one of the best tunneling spectra measured at T ~ A.SK 
in Bi-2212 [l^, as shown in the inset of Fig.l. 

Multi-gap model. - One can expect that the electronic 
inhomogeneity in HTSC may produce regions with a dis- 
tribution of gap amplitudes (A(z) and Ap{i)) and varia- 
tion in the local DOS. Recent STM and STS experiments 
on Bi-2212 and other high-Tc systems indicate that the 
gap inhomogeneities commonly exist in these materials 
regardless of doping level 22423. . Therefore, in order 
to reproduce the main features of the tunneling spectra 
of high- Tc cuprates, we have to consider the multi-gap 
case and the multi-channel tunneling processes, which 
contribute to the total tunneling current. At positive 
bias voltages V > 0, the tunneling of electrons from the 
normal metal into many regions of HTSC with different 
BCS DOS takes place and the resulting conductance is 



dl 



N- 



dV 



^ =^G4Al,(AT(^),ay) + ^2.(AT(^),ay)]. (7) 
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FIG. 3: Main panel: SIN tunneling spectrum measured on over- 
doped Bi-2212 at 43.1 K [27II fitted by using two-gap model, with 
A=31 and 18 meV; Ap= 22 and 15 meV. Left inset: fits of SIN 
tunneling spectra measured on underdoped Bi-2212 1271 by using 
three-gap model, with Ap=44, 28 and 20 meV and the set of gap 
values A=38, 26 and 17 meV for 46.4 K, A=37, 25 and 16 meV 
for 63.3 K and A=36, 24 and 15 meV for 76 K. Right inset: fit of 
SIN tunneling spectrum measured on slightly underdoped Bi-2212 
at 50 K 28] by using three gap model, with A=36, 24 and 15 meV; 
Ap= 73, 56 and 39 meV. 



In the case V < 0, the total current is the sum of tun- 
ncHng currents from many areas of HTSC with different 
local DOS {Dbcs{E, A{i)) and Df{E, Ap{i))) to the nor- 
mal metal. Then the resulting conductance is 

= J2 G,{AuiAT{i), -av) + v42,(At(*), -ay) 

i 

-f aF, (T) [Bh ( a; (i) , ) + ( a; W , ay )] } .(8) 

In such a multi-gap model, the tunneling spectra ex- 
hibit a more V-shaped behavior at low bias, the peak- 
dip-hump feature at negative bias and the asymmetry 
of the conductance peaks. With increasing temperature, 
the dip and peak on the negative bias side gradually dis- 
appear (see Fig. 2a), leaving the hump feature and the 
second conductance peak (on the positive bias side), as 
observed in tunneling experiments (l5j . Figure 2h shows 
that the conductance peaks become more asymmetric 
with increasing doping, as seen in experiments [13, . 

Comparison with the experiment. -The parameters en- 
tering into Eqs.© and ([5]) can be varied to fit experimen- 
tal data. The comparison of the theoretical results with 
the different experimental data on underdoped, slightly 
underdoped and overdoped Bi-2212 is presented in Fig. 3. 
We obtained the best fits to the experimental spectra by 
taking only two or three terms in Eqs.(l7]) and ([8]). In this 
way, we succeeded in fitting almost all of experimental 
conductance curves by taking different gap values. The 
V-shaped subgap feature, the asymmetric peaks and the 
dip-hump features, their temperature dependences ob- 
served in tunneling spectra of underdoped Bi-2212 (left 
inset in Fig. 3), slightly underdoped Bi-2212 (right inset 
in Fig. 3) and overdoped Bi-2212 (main panel in Fig. 3) 
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FIG. 4: Comparison of the tunneling conductance data on inhomo- 
geneous Bi-2212 (dashed line) 23) with the tunneling conductance 
calculated at 30 K (main panel) using the multi-gap model (A=52, 
45, 36, 27, 19 and 12 meV; Ap= 78, 65, 47, 34, 23 and 15 meV). 
Inset represents the conductance curve with pronounced shoulder- 
like feature, calculated at 30 K using the multi-gap model (A=59, 
43, 34, 25 and 16 meV; Ap= 64, 52, 44, 28 and 18 meV). 



are well reproduced. Moreover, the multi-gap model re- 
produces other tunneling spectra (with varying local gap 
value, ranging from 20 to 70 meV) and shoulders inside 
the conductance peaks observed in Bi-2212 [2^ In 
particular, this model reproduces rather well one of the 
experimental spectra of inhomogeneous Bi-2212 by 
taking six terms in Eqs.([7]) and ([5]), as shown in Fig. 4 
(main panel). Further, the conductance curve calculated 
using the multi-gap model (inset in Fig. 4) is similar to 
that in Fig.lc of Ref. [23| measured on inhomogeneous 
Bi-2212. 

We now discuss the relation between the BCS tunnel- 
ing gap and the SC order parameter. The unusually large 
reduced-gap values 2A/fcBTc ~ 7 — 22 observed in Bi- 
2212 [l^ compared to the BCS value 3.52 give evidence 
that the BCS gap determined by tunneling and ARPES 
measurements does not close at Tc and it is not related to 
the SC order parameter. While the peak suppression on 
the negative bias side near Tc observed in Bi-2212 is due 
to a spectral superposition of the tunneling conductances 
associated with the BCS DOS and quasi-free state DOS 
(originating from the polaron dissociation). The persis- 
tence of the conductance peak on the positive bias side 
well above Tc is evidence for the opening of a non-SC 
BCS gap at T* (for which the ratio 2/\/kBT* remains 
constant and close to the value 3.52 [Jll]). The pre- 
formed Cooper pairs condense into a SF Bose- liquid state 
at Tc (at which the SC order parameter appears) and the 
BCS pairing gap persists as the non-SC gap both below 
Tc and above Tc [l], 0] . The optical measurements (in- 
cluding tunneling spectroscopy and ARPES) are mainly 
sensitive to the excitation gaps at Ef, but such experi- 
mental probes compared with the thermodynamic meth- 
ods [13, HI] and the methods of critical magnetic field 
fiij and current [1^ measurements are insensitive to the 
identification of the SC order parameter as the SF con- 
densation energy or as the energy needed for destruction 
of a SF Bose-condensate (see also Refs. [H, 0, H^). 
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Conclusion - We have proposed a model describing the 
distinctive mechanisms of quasiparticle tunneUng across 
the SIN junction at negative and positive biases. The 
model incorporating effects of the BCS DOS and quasi- 
free state DOS (appearing at the polaron dissociation) 
at negative bias, and the gap inhomogeneity (i.e., multi- 
gap effects) reproduces the nearly U- and V-shaped and 
shoulder-like subgap features, peak-dip-hump structure 
and asymmetry of the conductance peaks and their evo- 
lution with temperature and doping as seen in tunneling 



spectra of Bi-2212. In this model, many unusual features 
of the timncling spectra observed in Bi-2212 on the nega- 
tive bias side arise from the spectral superposition of the 
tunneling conductances associated with the BCS DOS, 
quasi-free state DOS and multi-channel tunneling. 
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Baimatov, B.Y. Yavidov and B.V. Turimov for useful 
discussions. This work was supported by the Founda- 
tion of Uzbek Academy of Sciences, Grant No. FA-F2- 
F070-F075. 
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